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Helicases are molecular motors that unwind double-stranded nucleic acids (dsNA), such as DNA 
and RNA). Typically a helicase translocates along one of the NA single strands while unwinding 
and uses adenosine triphosphate (ATP) hydrolysis as an energy source. Here we model of a helicase 
motor that can switch between two states, which could represent two different points in the ATP 
hydrolysis cycle. Our model is an extension of the earlier Betterton-Jiilicher model of helicases 
to incorporate switching between two states. The main predictions of the model are the speed 
of unwinding of the dsNA and fluctuations around the average unwinding velocity. Motivated by 
a recent claim that the NS3 helicase of Hepatitis C virus follows a flashing ratchet mechanism, 
we have compared the experimental results for the NS3 helicase with a special limit of our model 
which corresponds to the flashing ratchet scenario. Our model accounts for one key feature of the 
experimental data on NS3 helicase. However, contradictory observations in experiments carried out 
under different conditions limit the ability to compare the model to experiments. 

PACS numbers: 87.16.Nn,82.39.-k,87.10.+c,87.15.Aa,05.40.-a,82.20.-w 



Helicases are enzymes that unwind double-stranded nucleic acids (dsNA) [l| . Helicase proteins typically translocate 
along one of the single strands and perform mechanical work while consuming chemical energy (usually supplied by 
the hydrolysis of ATP). Therefore, these NA translocases are molecular motors 0, Q which share common features 
with cytoskeletal molecular motors [1, [f| . 

All helicases undergo a biochemical cycle which typically involves ATP binding, ATP hydrolysis, and release of the 
hydrolysis products adenosine diphosphate (ADP) and and inorganic phosphate (Pj). An important question in the 
study of helicase mechanisms is to understand how the ATP hydrolysis cycle is coupled to the binding state and the 
motion of the helicase [g, |7fl. Helicases may exhibit changes in helicase/NA binding affinity when the helicase is bound 
to ATP, ADP/Pj, or neither; coordination of hydrolysis between different helicase subunits, and conformational 
changes in the helicase triggered by different steps in the hydrolysis cycle. Some helicases form hexamers (which 
include six ATPase domains), while others are members of the non-hexameric (dimeric or monomeric) group; different 
types of mechanochemical cycle have been suggested for the different structural classes d, Q. In all cases, one seeks 
to explain how the helicase coordinates NA binding and hydrolysis to move along single-stranded NA and unwind 
double-stranded NA. 

Here we develop a generic model of a helicase that switches between two biochemical states while translocating on 
ssNA. This is a simplified representation of the different states of the helicase during the ATP hydrolysis cycle. The 
model may be generally applicable to helicases for which the transition between two states is the key feature of the 
motion. In other words, this model should be a good approximation for helicases with more than two biochemical 
states if one transition is far slower than the others. We in corp orate such a two-state picture by extending the original 
Betterton-Jiilicher (BJ) model of NA helicases [Hf. 

Our work is also connected to two-state models that have been used extensively for a variety of molecular motors 
[l3l [lil . flEl . HH. Under a mean-field approximation, such models can be easily solved when periodic boundary 
conditions are imposed. However, the problem is usually more difficult with open boundary conditions. The model 
for helicase motion is even more complex because the position of one boundary (i.e., the ssNA-dsNA junction) varies 
randomly with time. Thus our work is also an extension of previous work on two-state models to the more difficult 
case of a fluctuating boundary. 

The two-state model developed here is consistent with the observation that binding and hydrolysis of ATP can 
modulate the affinity of a helicase for the nucleic-acid track [Tt], [H, EH • The flashing-ratchet mechanism suggested 
qualitatively for the hepatitis C virus non-structural protein 3 (HCV NS3) helicase [13, EH can be captured by a special 
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FIG. 1: Schematic of the model. The protein can exist in either of two chemical states (labeled 1 and 2) at each lattice site 
(labeled n). Sliding transitions (where n changes but the chemical state does not) occur at rate si/, etc., depending on the 
state and whether the transition is forward (toward increasing n) or backward (toward decreasing n). Chemical transitions 
(where the chemical state changes but n does not) occur at rates uj\2 (for the transition from 1 to 2) and u)2\ (for the transition 
from 2 to 1). Coupled transitions, where both the chemical state and n change, occur at rates 77 (for the transition from 2 
to 1 coupled to forward motion), (for the transition from 1 to 2 coupled to backward motion), q/ (for the transition from 1 
to 2 coupled to forward motion), and qb (for the transition from 2 to 1 coupled to backward motion). The nucleic acid single 
strand-double strand junction is at site m. The junction moves toward increasing m when the NA opens by one base (rate a) 
and toward decreasing m when the NA closes (rate j3). 

case of the generic model proposed here. In the flashing-ratchet [l3| picture, the motor protein switches between two 
states: one where the protein is tightly bound to the track, and another where the motor is weakly bound and can 
diffuse along the track. In this paper we make quantitative comparisons between our theoretical predictions for a 
passive helicase which follows the flashing-ratchet mechanism, and the experimental data for NS3 hclicasc. 

In section U we describe the ingredients of the model: the hclicasc, which can switch between two states and 
translocate on ssNA, and the fluctuating NA ss-ds junction. In section [TT| we calculate the single-strand translocation 
rate of the helicase. Section IIIII contains the model equations for double-strand unwinding, the transformation of the 
equations using midpoint and difference variables, and the general solutions for the velocity and diffusion coefficient. 
We describe the results for a hard-wall interaction between the helicase and junction in section HVl Using rate 
constants estimated from experiments on NS3 helicase, in section [V] we specialize to the flashing-ratchet scenario and 
make predictions specific to NS3. In section IVT1 we summarize our results. 

I. THE MODEL 

Here we develop a physical model for a helicase that moves on ssNA while cycling between two chemical states 
(labeled 1 and 2). Levin et al. suggested such a two-state model for NS3 helicase motion [13,[2l[. In this paper, we 
first consider a general two-state model, and later focus on the specific flashing-ratchet picture. 

In the traditional continuous models of Brownian ratchets, one first writes a Fokker-Planck equation. We use a 
discrete model, so our approach is based on master equations. The discrete approach can be useful when comparing 
to experiments. In the Fokker-Planck approach, one needs the explicit functional form of the fluctuating potential, 
which has not been measured for any real motor. In the discrete model, we bypass this difficulty by capturing the 
motor mechanism through a choice of rate constants (or transition probabilities) , many of which can be obtained from 
experiments (see section |V)) . 

In the discrete model, we represent the ssNA by a one-dimensional lattice where each site corresponds to a single 
base. We label each site by the integer index i. As in the BJ model @, we neglect the sequence inhomogeneity of the 
ssNA (in principle, the model can be extended to capture this feature, which may be important in some limits [22I]). 
The position of the helicase is denoted by the integer n. Most helicases have a fixed direction of translocation, either 
3' to 5' or 5' to 3' along the left-right asymmetric ssNA [g]. In our model the helicase translocates toward increasing 
n (from left to right in fig. [J). At any spatial position n, the helicase can be either in biochemical state 1 or 2. 

The model is fully described by the allowed transitions between states and the corresponding reaction rates. In 
general, we could have all transitions sketched in fig. [TJ Helicase "sliding" corresponds to transitions along the 
ssNA without a change in biochemical state of the protein. In state 1, these sliding transitions occur at rate sif 
(for increasing n) and sih (for decreasing n). When the helicase is in state 2, the forward/backward sliding rates 
are S2/ and S26- Physically, these transitions occur because of Brownian motion of the protein, decoupled from any 
biochemical state change. 

The helicase can undergo "chemical" transitions which correspond to a change in biochemical state without physical 
translocation along the ssNA. At fixed n, the rate of transition from state 1 to 2 occurs at rate 0J12 , while the reverse 
transition occurs at rate Finally, "coupled" mechanochemical transitions are those where a change of biochemical 
state and physical translocation occur together. If the helicase is located at n and is in state 2, then it can make a 
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transition to state 1 while moving forward to site n + 1 at rate r/; the corresponding reverse rate is r b . The transition 
of the helicase from state 1 to 2 while moving forward from n to n + 1 occurs at rate qj; the corresponding reverse 
rate is q b . 

If any of these reactions is coupled to ATP hydrolysis, then the forward/reverse transitions may be out of equilibrium 
and break the detailed balance relation. The Levin et al. model of HCV NS3 helicase suggests that ATP binding is 
required to remove the helicase from the tightly bound state [2fj, HH , implying that the 1 — > 2 transition at rate u>i 2 is 
determined by the ATP concentration. In the Levin et al. flashing-ratchet model, ATP hydrolysis and product release 
is coupled to the translocation and chemical transition back to state 1, which in our representation means that rates 
cu 2 i and Tf would be coupled to ATP hydrolysis and would therefore be out of equilibrium (see section |V| . 

The junction between ssNA and dsNA is labeled by m (see fig. [TJ. The dsNA opens and closes due to thermal 
fluctuations. When the helicase and junction are far apart, the opening rate is a and the closing rate 0. We assume 
that these rates are independent of the NA base sequence and that the only fluctuations are those for which the NA 
opens or closes at the ss-ds fork. Following the BJ model @, we neglect the possibility of any jump > 1 bp in the 
position of the ssNA-dsNA junction. However, this approximation is justified because, at the temperatures of our 
interest (i.e., sufficiently below the melting temperature of the dsDNA) the spontaneous formation of bubbles is rare. 
Since the NA breathing results from thermal fluctuations, the rates a and (3 satisfy detailed balance: % = e~ AG , 
where AG is the free energy of one base-pair bond in units of kT. 

The main quantity of interest is the speed of unwinding of dsNA by a helicase. We derive an analytical expression 
for the unwinding velocity. We compare the predicted velocity with the corresponding experimental data for a specific 
helicase, NS3 helicase of hepatitis C virus. Although we also derive an analytical expression for the diffusion constant 
of the helicase, we do not compare it with experimental data for any specific helicase. 

In this work we analyze passive unwinding, which is equivalent to a hard-wall interaction potential in the BJ model 
[To| . In passive unwinding, the helicase acts as a block to NA closing when adjacent to the junction. The protein 
moves forward only when thermal fluctuations open a basepair at the NA ss-ds junction. This means that when the 
helicase and junction are adjacent (j — 1), the helicase cannot hop forward (all helicase forward rates, s\f(j = 1), 
s 2f{j = 1), Tf(j = 1), and qf(j = 1), are zero) and the NA cannot close = 1) = 0). Otherwise, the rates are 
unaffected by the helicase-junction interaction. 



II. SINGLE-STRAND TRANSLOCATION 



In order to motivate our approach, we first formulate the equations for a helicase sufficiently far from the ssNA-dsNA 
junction so that it translocates on ssNA without any dsNA unwinding activity. Let Vp(n,t) denote the probability 
that, at time t, the helicase is located at site n and is in the chemical state [i. We will drop the reference to the time 
dependence of V^in). The master equations governing the time evolution of V^in) are 



rfPi(n) 
dt 

and 

dP 2 (n) 



-(wi2 + si/ + si 6 + 9/ + n)Vi{n) + si f P 1 (n - 1) + r f T> 2 {n - 1) 

s lb V 1 {n + 1) + q b V 2 {n + 1) + io 21 V 2 (n), (1) 



dt 



-(w 2 i + s 2f + s 2b + Tf + qb)V 2 {n) + s 2f V 2 (n - 1) + q f Vi{n - 1) 
+ s 2b V 2 {n + l)+r b Vi{n+l)+u) l2 Vi{n). (2) 

Summing these equations, we find the total probability V(n) — Vi(n) + V 2 (n) satisfies 

dV{n) 



dt 



= ~( s if + s u + qj + n)Vi{n) - (s 2 / + s 2b + r f + q b )V 2 {n) + (si/ + qf)V\{n - 1) + (s 2 / + r f )V 2 {n - 1) 
+ (si b + r b )Vi(n+l) + (s2b + qb)V 2 (n+l). (3) 



These equations have a translationally invariant steady-state solution where V^n) is independent of n. In this case, 
we expect that the probability in state 2 is a multiple of the probability in state 1: 

Pa(n) = aVi{n), (4) 

which means that V[ri) = (1 + a)Vi(n). 
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In this case, the master equation for the total probability can be written as a hopping model with effective rates 
kf for forward transitions and k b for backward transitions. At steady state, 



where 



and the expression 



= kfV(n-l)-(kf + k b )V(n) + k b V(n+l), (5) 



_ si/ + q f +cr(s 2 f +Tf) , . 

k f - ^ ' [ > 

, 8 a + r b + a(s 2b + qb) frTS 

h = — . (7) 



^12 + g f + n 

r f +Qb + ^21 ' 



has been obtained from cqn. ([T]) at steady state, assuming translational invariance. The mean single-strand translo- 
cation velocity is v ss — kf — kb- 

III. DOUBLE-STRAND UNWINDING: MODEL EQUATIONS 

In this section we extend the formulation of the preceding section by incorporating helicase-catalyzed dsNA un- 
winding. Let Vfj.(n,m;t) denote the probability that, at time t, the helicase is at located at n and is in the chemical 
state fi, while the ss-ds junction is at to. We will drop the reference to the time dependence of Vfi{n, m). The master 
equations governing the time evolution of V^in, to) are given by 



dVi{n, m) 
Jt 

and 

dV-iin, m) 



-{a + (3 + W12 + sif + su + qf + r b )Vi(n,m) + sifVi(n — l,m) + rfV^in — l,m) 

s\bPi{n + l,m) + qbP-xin + 1, m) + UziP^n, m) + aVi(n, m — 1) + /3Pi(n, m+ 1) (m > n). (9) 

-(a + (3 + lu 2 i + s 2 f + s 2 b + r f + qb)V 2 {n,m) + s 2f V 2 (n - 1, m) + qfVi(n - 1, m) 

s 2 bP 2 {n + 1, to) + r b Vi(n + 1, m) + 07i 2 'Pi(n, m) + aV 2 (n, to — 1) + (3V 2 (n, m + 1) (m > rt). (10) 



Note that the rates depend on the separation m — n; this notation is omitted for clarity. We assume the interaction 
potential is the same for both chemical states, so that the position-dependent NA opening and closing rates a and (3 
are independent of the chemical state. 

Next we change variables to work with the difference j = m — n and midpoint I = 21 = m + n positions of the 
helicase-j unction complex. Rewriting eqns. © and (|10p we have 



dt 



-{a + [3 + wia + si f + s u + q f + r b )Vi(j, I) + s lf Vi{j + 1, 1 - 1) + r f T 2 {j + 1,1-1) 
sibPi{j -1,1 + 1) + q b V 2 {j -1,1 + 1)+ uj 2 iV 2 (j, I) + aPi{j -1,1-1) + (3Vi{j + 1,1 + 1) 

(J > 0). (11) 

-{a + 13 + oj 21 + s 2f + s 2b + r f + q b )T 2 (j, I) + s 2f P 2 (j + 1,1-1) + q f Vi(j + 1,1-1) 

+ s 2b V 2 {j -1,1 + 1) + r b Vi(j -1,1 + 1) + uJi 2 Vi(j, I) + aV 2 {j -1,1-1) + /3V 2 {j + 1,1 + 1) 

(J > 0). (12) 

Again, the rates vary with j. However, the rates are independent of /, so we can sum over the position of the complex 
center of mass: 



and 

dP 2 U,l) 
dt 



iH?') = £W) (is) 
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Applying the sum over I to eans. fTTj ) and (fT2|) we find 



dt 



{a + p + Lu 12 + si f + s lb + q f + r 6 )Pi(j) + {si f + /3)Pi(j + 1) + r f P 2 {j + 1) 



+ (si6 + a)Pi(i-l) + g 6 P2(i-l)+W2iP2(j)- (14) 

and 



dt 



= -(a + /3 + w 2 i + s 2/ + s 2b + r/ + q b )P 2 (j) + (« a / + /9) p 2(j + 1) + q f Pi(j + 1) 



+ (s 2b + a)P 2 (j-l)+r b Pi(j-l)+u 12 Pi(j). (15) 
We consider the total probability by summing eqns. (|14| and (|15[) . Defining the total probability current 

= ~ W + 1) + + ^)^i(j) + (S2b + q b )P2(j) - (sif + q f )Pi(j + 1) - (s 2f + r f )P 2 (j + 1), (16) 

the total probability satisfies 

dP(j) 



dt 



= -I(j)+I(j-l). (17) 



At steady state P(j) is time independent, so = I(j — 1). Further, since U(j) — » oo as j — > — oo, this constant 
probability flux must be zero, i.e., — for all j. 

Adding the two eqns. (fTTj) and (TH?)) and defining V(j,l) = Vi(j,l) +V 2 (j,l), we get 

= -(a + /3)P(j, + aV(j - 1, / - 1) + /3P(j + 1, 1 + 1) + (s lf + q^V^j + 1,1-1) 
+ (r f + s 2f )V 2 {j + 1,1-1) + (s lb + r b )Vi{j -l,l+l) + (q b + s 2b )V 2 (j - 1, 1 + 1) 

+ oj 2 iV 2 (j, l) + LJuViij, l) - (^12 + si/ + si 6 + qf + n)V\{j, l) 

- (uJ 2 i+s 2f + s 2b + r f +q b )V 2 (j,l). (18) 
The probability distribution in I at time t is 

n(Z;i) = £>(;M;i) (19) 

3 

Note that, by definition, U(l;t) is independent of the chemical state of the helicase For times much longer than the 
relaxation time of the difference variable j, we can assume 

V tl (j,l)=P f ,(j)U(l) ( M = lor2) (20) 

Starting from the eqn. (|18p . one can derive 

= yji(i _ i) _ ( u + W )U(1) + wll(l + 1) (21) 

where 

u = J2 aP(j) + (s lf + q f )Pi(j) + (s 2f + r f )P 2 (j), (22) 

3 

and 

w = PP(j) + (*16 + n)Pl(i) + (s 2b + q b )P 2 {j). (23) 

Thus the motion of the helicase-j unction complex is a combination of drift and diffusion. Note that in the special case 
u = w the drift vanishes and the dynamics of I becomes purely diffusive. 
As in ref. [l(J, the average speed of unwinding is v — h (u — w) , or 

v = \ Y.( a ~ ® P W + ( Sl / + 9/ - Slb _ r b) p iU) + (s2 f +r f - s 2b - q b )P 2 (j). (24) 

3 

Similarly, the diffusion coefficient is D = j(u + w), which is 

D = \ + + ( Sl / + q f + Slb + rb ) p i(-?) + (*»/ + r / + s 2b + ?fe)P2(i). (25) 

3 

Note that if the sliding transitions represent unbiased diffusion, then the forward and backward rates s^f and s^ b are 
equal. Then the terms involving the sliding rates drop out from the expression for v but not from that for D. 
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IV. SOLUTION 

In order to evaluate the expressions for the unwinding velocity and diffusion coefficient, we must determine P\{j) 
and P 2 (j) hi terms of the rate constants. Consider the result of summing eqns. (|14[) and (|15p over j to determine 
equations for the total probability of being in state 1, Pi and the total probability of being in state 2, P 2 . We can 
write these equations as 

= -*>uPi + *to. p *> (26) 

at 

= -k 21 P 2 + k 12 P u (27) 

where the rate constant ki2 depends on q/, and r b , and &21 depends on W21, »"/, and q b . The steady-state solution 
has P 2 = k 12 /k 2 iPi. 

This observation suggests a translationally invariant solution for P\(j) and P 2 {j) when the rates are constant. We 
consider the case where the relative probability of being in state 1 or 2 is translationally invariant (independent of j). 
This must occur if the hopping rates are constant or spatially vary in the same way (for example, if states 1 and 2 have 
the same interaction potential with the dsNA). Since we are primarily interested in a passive helicase with constant 
rates, we will focus on this case. Because of the translational invariance, the probability in state 2 is a multiple of the 
probability in state 1, so that 

Pi(J)=lPi<J)> (28) 
The zero-current relation requires that eqn. (|16p equal zero, which requires 

(/J + si f + q f )Pi(j + 1) + (fi + s 2f + r f )P 2 (j + 1) = {a + s lb + r 6 )Pi(j) + (a + s 2b + q b )P 2 (j). (29) 

We can plug in to eqn. (|2"5|) and solve for the unknown constant 7. We can rewrite eqn. (|2"9"1) as a recursion relation 
that relates Pi(j + 1) to P\{j): 

PiU + l) = «(1 + 7) + si6 + r 6 + 7(326 + 96) = c (m) 
PiU) 0(1 + 7) + sif + q f + 7(s 2/ + r f ) 

Note that c is a function of 7. While it is possible to solve coupled equations for c and 7 in general, the resulting 
expressions are long and not useful for developing intuition. Instead, we use the approximation relevant for helicases 
that a and f3, the opening and closing rates of the NA, are several orders of magnitude larger than the other rates 
in the problem (see reference [13] , where experimental data from reference [23| was used to estimate the opening rate 
a ~ 10 7 s _1 ; other rates in the problem are of order 10 2 s" 1 ). In this case, eqn. ([30]) reduces to 

c«f- (3D 

Throughout the remainder of this paper, we will use this approximate value of c. Note that because a and f3 are 
constant, c is also constant and eqn. l[30|) shows that Pi(j) has power-law decay with increasing j (as in the BJ model 
for a passive helicase [Toj). 

Using eqns. (|28|) and (|30|) in eqn. (T14]) at steady state, and imposing the requirement that P\(j) cannot vanish for 
arbitrary j, we find a unique expression for 7: 

s i/(l -°)- SibjcT^ -l)+r b + q f + u 12 , qoA 

7 = : — — : • (32) 

erf + c 1 q b + oj 2 i 

With this result, we can evaluate eqns. and ([23]) and express v and I? in a fashion analogous to the expressions 
in the simpler BJ model: 
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if we define the effective rates 

a = a(l + 7), (35) 

b = 0(1+7), (36) 

k+ = 7 (s 2/ + r/ ) + si/ + <?/, (37) 

= 7( s 26 + 3b) + si6 + n>. (38) 

Next we evaluate the sums in eqns. (|33p and lj3~i)) , noting that P\{j) = P\c? and taking into account that for j = 1 
the rates fc + and 6 are zero. The result is 

ck~^~ k . , . 

u = WW (39) 

a ck + + k~ , . 

D= 2 + IoTtT (40) 

Equations (|39|) and (J^Dj) are the main results. 

Note that under most conditions the NA opening and closing rate a is orders of magnitude larger than the other 
rates, and therefore D w a/2. 



V. COMPARISON WITH NS3 HELICASE 



The NS3 helicase of the hepatitis C virus (HCV) is important for HCV replication, and is therefore a potential 
drug target [HJ|. NS3 is also an interesting model helicase because it is the only currently known helicase capable of 
unwinding both dsRNA and dsDNA (2f| Hy] . The flashing-ratchet mechanism proposed for NS3 helicase in ref . 21] is 
a special case of the two-state model which we have developed in the preceding sections. In this section, we first briefly 
summarize the experimental data on NS3 helicase and their mutually contradictory interpretations which highlight the 
current debates in the literature. Then, we present analytical results for the special case of our model which captures 
the flashing ratchet mechanism. We compare these theoretical predictions with the corresponding experimental data 
for NS3 helicase. The comparisons are, however, limited by the contradictions between the observations in different 
experiments, many of which have been performed under different conditions. 



A. Summary of experimental results on NS3 helicase 



To compare our model to experiments on NS3 helicase, we would ideally like to know the enzyme step size, the 
single-strand translocation rate, and the double-strand unwinding rate — including information on how it varies with 
NA sequence or applied force. Interpretation of experimental data on NS3 is complicated by differences in experiments 
done by different research groups. Some groups study the full-length NS3 protein, including the helicase and protease 
domains [H HI H! HS HI, while others study the helicase domain only [H HH HI, S3, HI, H HI]. Moreover, 
genetically different versions of NS3 can have different properties (36| . The NS3 protein can also function in different 
oligomeric states. In bulk solution experiments, full-length NS3 seems to function best as a dimer or higher-order 
oligomer (3?| , but single- molecule experiments can observe unwinding by NS3 monomers [13, HH . The helicase domain 
NS3h appears not to form dimers in solution [TH, H3, H|[, but multiple copies of the protein can bind to ssNA and 
unwind dsNA [33]. In at least one experiment, the kinetic parameters did not vary with the length of the ss tail used 
to load NS3h, suggesting that the helicase mechanism may not depend on whether the protein is a monomer or dimer 

M 

Contradictory claims have been made in the literature on the qualitative description of NS3 helicase as well as 
on its quantitative characteristics. First, we consider the empirical evidence for the stepping pattern and the step 
size of NS3 helicase. Recently a detailed computational model of NS3, based on known crystal structures, supported 
the idea of single-base "inchworm" motion taken by NS3 monomers. This model of Zheng et al. proposes a major 
protein conformational change which is triggered by ATP binding and is coupled to forward motion of the helicase 
[39l ] . Models based on structural studies of NS3 have suggested single-base steps [H, H(| ■ Similarly, structures of the 
distantly related Hel308 helicase, which shows some structural similarities to NS3, supports the idea of a ratchet-like 
mechanism during the ATP cycle [4lj . However, most experimental efforts to determine the step size don't support 
single-base steps. Bulk kinetic experiments have given a kinetic step size of 9-17 basepairs, depending on protein form 
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FIG. 2: Schematic of the simplified model that represents a flashing ratchet. The protein can exist in either of two chemical 
states (labeled 1 and 2) at each lattice site (labeled n). Sliding transitions (where n changes but the state does not) occur 
only in state 2 at rate S2- Chemical transitions (where the state changes but n does not) occur at rates W12 (for the transition 
from 1 to 2) and U21 (for the transition from 2 to 1). A coupled transition (where both the state and n change) occurs at rate 
77 (for the transition from 2 to 1 coupled to forward motion). The nucleic acid single strand-double strand junction is at site 
m. The junction moves toward increasing m when the NA opens by one base (rate a) and toward decreasing m when the NA 
closes (rate f3). 



and unwinding substrate [U H^, [35| . Single- molecule experiments on monomers of full-length NS3 have suggested 
a step size of 11 basepairs with 3 basepair substeps 30] or 3 basepairs with 1 basepair substeps [31| . The most 
recent single-molecule work has proposed that the fundamental step size is one basepair, with pauses occurring less 
frequently as part of the ssNA bound to the helicase occasionally "rips" off [3l[ . 

Next we summarize the current estimates of ss translocation rate and the speed of double-strand unwinding by NS3 
helicase. The maximum ss translocation rate can be estimated from experiments that measure the ATP hydrolysis 
rate. In one experiment, the NS3h rate of ATP hydrolysis had a maximum k cat of 80 s" 1 in the presence of the 
single-stranded oligo dUis (Hf. Assuming that during ss translocation the helicase hydrolyzes 1 ATP per step, this 
measurement sets an upper bound on the ss translocation velocity of 80 bases s _1 . The double-strand unwinding 
velocity of NS3 has been estimated from bulk and single-molecule experiments. In one single-turnover bulk kinetic 
study, the maximum unwinding rate of NS3h was 2.7 bp s _1 (3||; similar results were found by another group [28| . 
Full-length NS3 may unwind at higher velocities, up to 16.5 bp s _1 [1?], HI]. In single-molecule experiments with 
applied force, full-length NS3 monomers unwind at force-independent rates of 50 bp s _1 [301] . This relatively high 
velocity may be possible because of the applied force that reduces the energetic cost of opening the NA. In single- 
molecule FRET experiments on full-length NS3 monomers where no force is applied, an unwinding rate of k « 0.9 
s _1 was measured for one base pair substeps [3l[ — a value closer to the bulk value measured for NS3h. 

Finally, we examine the experimental data to investigate whether the unwinding by NS3 helicase is active or 
passive. The dependence of the unwinding rate on the base-pair binding free energy was measured both in single- 
molecule and bulk experiments. In the work of Dumont et al., the RNA unwinding rate of full-length NS3 monomers 
was approximately independent of applied force in the range 9-17 pN [3fjJ. In this experiment, the applied force 
was relatively high: the double strand melted at a force of 20 pN. In single-molecule experiments using a similar 
experimental setup, Cheng et al. [HJ observed a significant effect of varying the RNA sequence on the NS3 unwinding 
rate. This observation of Cheng et al. indicates that a passive unwinding mechanism may not be adequate to explain 
the behavior of full-length NS3 helicase. Further, the apparent contradiction between the observations of Cheng et 
al. [42| and Dumont et al. [30j may be reconciled if we abandon the simple physical picture in which the base-pair 
binding free energy can be altered in a similar way by applied force or by changing the sequence. Recent bulk 
measurements examined the effects of sequence variation on the unwinding rate of NS3h [431 ] ; this work is discussed 
below where we compare our theoretical predictions to experimental results. 

In order to motivate our minimal model for the NS3 helicase, we now discuss the affinity of NS3h to NA and its 
modulation during the ATP hydrolysis cycle. Binding experiments on NS3h found that when the helicase is bound to 
an ATP analogue, it binds to NA more weakly than when not bound to ATP or ADP [33|, HH ■ The change in binding 
free energy is approximately 6 kT at room temperature (15 kJ mol -1 ) [2l|. In addition, the affinity of NS3h for ADP 
is low, so release of hydrolysis products is expected to be rapid [3IJ . These observations are the basis of the proposed 
flashing-ratchet mechanism of_NS3h. (However, we note that another work has found no dependence of NA binding 
on the ATP hydrolysis state [28(; the source of this difference is unclear.) 



B. Flashing-ratchet model of NS3 helicase 

Here we consider a special case of our model which corresponds to a flashing ratchet mechanism. Levin et al. 
posed that NS3 helicase switches between two states: one tightly bound to the ssNA, the other weakly bound [13, 
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This scenario is referred to in the physics literature as a flashing ratchet [HI]. When applying the flashing ratchet 
scenario to NS3, the tightly bound state is represented by a periodic sawtooth potential (with periodicity of one 
ssNA base pair) and the weakly bound state is represented by a uniform (weakly position- independent) potential [2l[ . 
When comparing to the flashing-ratchet scenario, we will consider state 1 to represent the strongly bound (S) state 
and state 2 the weakly bound (W) state. By comparing the theoretical predictions for this special case of our model 
with the experimental data for NS3 helicase, we test whether or not NS3 follows the flashing ratchet mechanism. 

We assume that no sliding is possible in the tightly-bound state 1, so s±f = sib = 0, and that the sliding is unbiased 
in state 2, so S2f = S2b = *2- To connect with the flashing-ratchet scenario and for simplicity, we assume that the 
rates qf = qb = rh — (see fig. [2j . With these assumptions, we find that the rate of ss translocation is (from eqns. (J6]) 
and 0) 

V ss = Wi 2 , f , , (41) 

Tj + Wl2 + W21 

and the rate of ds unwinding is 

^12 {erf - (1 - c)s 2 ) 



v u = 



2 erf + LU 12 + L021 



(42) 



The excitation rate u>i2 is associated with ATP binding, and so is assumed proportional to ATP concentration. 
Therefore we write W12 — w D [ATP]. The rates W21 and r/ represent the relaxation from the weakly bound to the 
tightly bound state that occurs after ATP hydrolysis, product release, and diffusion in the weakly bound state. For 
a flashing ratchet, a high rate of forward motion will occur when the positions of the energy barriers and the time 
constants are such that forward movement (rate 77) and return to the same place after one cycle (rate UJ21) occur 
with equal probability. To match this optimal case, we therefore assume that U21 = r/. Further, we assume that the 
sliding rate S2 is small compared to the other rates; for concreteness we will suppose S2 = erf with e = 0.1 unless 
otherwise stated. The velocities then become 



77^ [ATP] 
w [ATP] +2r/' 
(c-e(l-c)) r/i^ATP] 

2 w [ATP] + (l + c)r/' 



(43) 
(44) 



Both v ss and v u are consistent with the Michaelis-Menten equation for enzyme kinetics, but with slightly different 
forms. Their ratio is 

= (c-e(l-c)) q; [ATP] + 2r f 
Vss 2 w [ATP] + (l + c)r/ 1 ' 

In other words, we predict that the ratio of the unwinding velocity to the single-strand translocation velocity depends 
on ATP concentration. If we average over sequence variation in DNA Q, we get the estimate c = a/0 ~ 1/7. For the 
purpose of quantitative illustration of the variation of — with ATP concentration, let us assume e = 1/10. Then, 
2*- w 0.029 at high ATP concentration and w 0.05 at low ATP concentration. This suggests that the ratio of the 
unwinding velocity to the single-strand translocation velocity could vary significantly with ATP concentration — the 
change is almost a factor of 2 for this example. 

Next, we estimate v u and v ss for NS3 helicase. The single-strand translocation and unwinding velocities are fully 
determined by the parameters c, rj, uj 0i e, and ATP concentration; we now extract estimates of r/ and ix> from 
experimental data. In experiments at high ATP concentration and in the presence of ssNA, NS3h shows a maximum 
ATP hydrolysis rate of 80 s _1 32]. If we take this value as the limiting ss-translocation rate and assume single 
base-pair steps, then v ss — 80 nt s _1 in the limit of high ATP concentration. Using this estimate of v ss in eqn. (I43[) . 
we get the estimate rf = 80 s _1 . This, in turn, implies that at high ATP concentration the unwinding velocity 
v u ~ 0.029w ss w 2.3 bp s _1 . This value is comparable to the values of 2.7 bp s _1 f35| found for NS3h and 0.9 bp 
s _1 found for the one-bp substeps of full-length NS3 [31]. We note that the unwinding velocity v u <C v ss , as should 
be expected for this model which assumes a passive helicase mechanism. Experiments studying how NS3 ATPase 
activity [32| and unwinding (30j vary with ATP concentration found a similar Michaelis constant K m w 90 /j,M. Using 
this value of K m in eqn. (|4*3"|) , we estimate uo — 2rf/K m « 1.8 fiM^ 1 s _1 . 

The only remaining unknown parameter is e = S2/?"/, ^ ne ratio of the sliding rate to the forward transition rate. 
A smaller value of e means that the sliding transitions in the weakly bound state are less probable (see fig. [5]) . A 
higher value of e means that sliding transitions in the weakly bound state are more probable. This parameter has an 
important effect on the dependence of the helicase velocity on the base-pair binding free energy. 
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FIG. 3: Dependence of the unwinding velocity on the base-pair binding free energy. The reference state is a value c = 1/7, 
which represents a sequence- averaged value for DNA. The additional destabilization energy AG (in units of kT) represents a 
free energy change that favors NA opening. When e increases, the dependence of the velocity on AG becomes more pronounced. 
However, decreasing e can not flatten the curve indefinitely. 



To study the effects of varying the base-pair binding free energy, we focus on the limit of high ATP concentration. 
In this case, if AG is the free energy of destabilization of base-pair binding, the parameter c = a/ (3 varies according 
to c = c a e AG . Therefore, at high ATP concentration, the unwinding velocity varies as 

lim v u = r s . (46) 

[ATP]^oo 2 1 

The unwinding velocity increases exponentially if the NA is destabilized, as one would expect for a passive helicase. 
However, the precise shape of the curve of unwinding velocity versus AG depends on e. In the limit e — > 0, which 
physically means no helicase sliding transitions occur in the weakly bound state, the unwinding velocity varies with 
AG as a simple exponential: 

lim Vu = c ^LL e ^_ (47) 

As e increases, the helicase can slide in the weakly bound state. This allows more rapid unwinding by the helicase: 
when the dsNA is destabilized, the ds base just ahead of the helicase has an increased probability to be open. Rather 
than wait for the helicase chemical transitions to move forward, the helicase can take advantage of this increased 
junction open probability and slide forward. This allows the steeper rate of increase of v u with AG seen in fig.[3l This 
prediction is qualitatively consistent with the result of Tackett et al. who found that full-length NS3 unwound 
double strands with higher melting temperatures less efficiently. However, in the single-molecule experiments of 
Dumont et al. the unwinding rate of full-length NS3 monomers was practically independent of applied force in the 
range 9-17 pN [3(3]. This disagrees with the prediction of this model, if the only effect of the applied force is to 
change the binding free energy per base pair. However, this physical interpretation is clearly not valid, because recent 
experiments from the same lab find a significant variation in the RNA unwinding rate of full-length NS3 with the 
variation of the base composition of the RNA [12] • Reconciliation of the apparent contradictions in these experimental 
observations is possible by assuming an active helicase mechanism which, however, is not incorporated in the current 
version of our model. Analyzing data from bulk experiments, Donmez et al. (43j claimed that the the variation of NS3h 
unwinding velocity with base-pair binding free energy is inconsistent with a passive helicase mechanism. However, 
this conclusion is drawn from an analysis based on a reported single-strand translocation velocity of 6.4 bases s _1 , 
which is much lower than the value of 80 bases s _1 mentioned above. A ss translocation rate of 80 bases s _1 is 
an upper limit, assuming the helicase hydrolyzes 1 ATP per single-base step. If the helicase on average hydrolyzes 
> 1 ATP per step, the ss translocation rate would be lower. A lower ss translocation rate would lead to an even 
larger disagreement between the passive helicase model we presented and the experimental data. We believe that a 
conclusive comparison between our model of a flashing-ratchet mechanism for NS3 helicase and the experimental data 
is not possible because of the contradictory reports of experimental studies. 
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VI. CONCLUSION 

In this paper we have developed a general model of unwinding of a double-stranded nucleic acid molecule by a 
helicase motor. To capture some of the key features of the helicase mechanochemical cycle, we have modeled helicase 
switching between two chemical states. In this model, the sites of a discrete lattice represent the positions of the 
individual bases on the ssNA. At any spatial position, the helicase can exist in either of the two allowed chemical 
states. This model should be generally applicable to helicases where one of the transitions in the mechanochemical 
cycle is much slower than the other transitions. In this work, we have considered only a passive helicase mechanism — 
the helicase at the junction must wait for thermal fluctuations to open the dsNA before it can advance. In future 
work, it would be valuable to extend the model to include active destabilization of the dsNA by the helicase. 

To compare the model in detail to experimental data, we focused on a special case which captures the flashing- 
ratchet mechanism proposed for the NS3 helicase [2l|. Solving the master equations for this model at steady state, 
we have calculated the speed of unwinding and the speed of single-strand translocation. The ratio of the unwinding 
velocity to the ss translocation velocity varies with ATP concentration as well as with the base-pair binding free 
energy. 

Our comparison to experimental data on NS3 helicase suggests that the model captures some features of the 
experiments. However, the experimental literature on NS3 contains contradictory results. This may be a result of the 
different genetic variants, protein truncations, oligomeric states, substrates, and buffer conditions used by different 
laboratories. A set of detailed experiments by different labs under consistent conditions may be important to fully 
understand the unwinding mechanism of NS3 helicase. 
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